Dissipation of Quantum Turbulence in the Zero Temperature Limit 
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Turbulence, produced by an impulsive spin-down from angular velocity SI to rest of a cube-shaped 
container, is investigated in superfluid 4 He at temperatures 0.08 K - 1.6 K. The density of quantized 
vortex lines L is measured by scattering negative ions. Homogeneous turbulence develops after time 
t w 20/J7 and decays as L oc t~ 3 ^ 2 . The corresponding energy flux e = v'(kL) 2 oc t~ 3 is characteristic 
of quasi-classical turbulence at high Re with a saturated energy-containing length. The effective 
kinematic viscosity in T — limit is v' = 0.003k;, where k — 10 



3 cm 2 s 1 is the circulation quantum. 



PACS numbers: 67.40.Vs, 47.27.Gs, 47.32.-y, 47.37.+q 

We report the first measurements of the decay of ho- 
mogeneous turbulence in superfluid 4 He in the zero tem- 
perature limit, where thermal excitations are effectively 
absent. The turbulence, produced by bringing a cubi- 
cal container rapidly to rest from a state of rotational 
equilibrium, has the classical Kolmogorov form on large 
length scales, even though on short scales rotational mo- 
tion is severely restricted by the quantization of circula- 
tion. At T < 0.8 K, we observe a drop in the effective 
kinematic viscosity, by a factor of ~ 30, to an ultimate 
T = value of 0.003k. We associate this drop with a 
regime in which the energy flux in the Kolmogorov cas- 
cade, when it reaches the quantum scale, must be trans- 
ferred to a Kelvin-wave cascade in which energy is car- 
ried to much smaller, dissipative, length scales by Kelvin 
waves on the individual quantized vortex lines. 

Superfluid He (strictly the superfluid component) is 
inviscid, and rotational flow can be achieved only with 
line defects, quantized vortices each carrying one 
quantum of irrotational circulation k = h/m outside a 
core of radius ao ~ 1 A. Turbulence in the superfluid 
( "quantum turbulence" ) can exist as a dynamic tangle of 
these vortices [2] , characterized by their total length per 
unit volume, L. On length scales much larger than the 
average distance between vortices, I = L -1 / 2 , rotational 
flow is possible as the result of partial alignment of the 
vortices [!, H| . Hence, at sufficiently large scales the su- 
perfluid should be able to support quasi-classical turbu- 
lence, similar to that in classical liquids at large Reynolds 
numbers, in which the energy cascades towards smaller 
eddies until it is dissipated [5|. Numerical simulations 
:6j confirm this. Dissipation is through the interaction 
of quantized vortices with excitations in the liquid. In 
the zero-temperature limit emission of phonons by Kelvin 
waves at wavelengths ~ 100 A is expected [1, Q, while 
above T as 0.5 K, scattering of the thermal excitations 
("mutual friction") should dominate [3J]. 

We see that in the T = limit there is a need to 
transfer energy at scales of order t from quasi-classical 
turbulent flow to ultra-quantum Kelvin- wave turbulence. 
L'vov et al. Q have suggested that this transfer involves 
a bottleneck, so that continuity of the energy flux requires 



that the vortex line density at mesoscales ~ £ should be 
enhanced. Kozik and Svistunov Q argue that reconnec- 
tions operating in a range of wave numbers, covering an 
intermediate regime between the Kolmogorov spectrum 
and the Kelvin- wave spectrum, allow the bottleneck to 
be effectively bypassed. 

The rate of dissipation, mediated by vortex lines, per 
unit mass can be written in the form 

e = v\kL)\ (1) 

where n 2 L 2 is an effective total mean square vorticity [ij 
and the parameter v' is an "effective kinematic viscosity" . 

In previous studies of quasi-classical turbulence in su- 
perfluid 4 He at T > IK, turbulence that was nearly 
isotropic and homogeneous was gene rated by either ro- 
tating blades [l(| or a towed grid [111]. To extract v 1 [ill ], 
one can monitor the late-time decay of the vortex density, 
which has the form L(t) oc v'~ 

-i/2 t -3/2^ xhe f ii 0W i nK 

assumptions lead to this type of decay. The turbulent 
energy should be concentrated in the largest eddies char- 
acterized by a size d and velocity u, so that the total tur- 
bulent energy per unit mass is E « u 2 /2; d must be con- 
stant in time, as is the case if it is limited by the container 
size; and the lifetime of the largest eddies with respect to 
breaking up into smaller eddies without dissipation is of 
order the turn-over time r ~ d/u. Hence, the energy flux 

e = -E ~ u 3 /d E 3 / 2 /d. This yields E ~ d 2 (t + t*y 2 

and e ~ d 2 (t + t*)~ 3 , where t = coincides with the time 
of activation. Here, t* is an arbitrary constant absorb- 
ing the duration of the transient processes and the initial 
value of the turbulent energy. It is usually of order r, and 
hence in what follows we will omit it as we are interested 
in the late-time decay for i > t. Finally, using Eq. (1) 
we obtain kL ~ dz/ -1 / 2 ^ 3 / 2 . We shall use this relation 
with the numerical prefactor, 



L= (3C) 3 / 2 (2^ K )- 1 ^'- 1 / 2 <- 3 / 2 , 



(2) 



derived by Stalp et al. [ll|, [T2[ for homogeneous isotropic 
turbulence in a channel of square (d x d) cross-section, 
using the Kolmogorov energy spectrum with C — 1.5, 
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and with an assumed sharp cut-off for k < 2w/d. 

Bradley et al. [l3[ have measured L(t) produced by 
a vibrating grid in 3 He-B down to T = limit, and ob- 
served a free decay L(t) oc t~ 3 ^ 2 . By assuming a value of 
d in Eq. (1) equal to the measured spatial spread of the 
turbulence, they inferred v' ~ 0.2k. While their obser- 
vations are interesting and suggestive, the turbulence is 
not homogeneous, and the size of the energy-containing 
eddies may differ from the spatial extent of the turbu- 
lence, so that the value of v' is ambiguous. Other exper- 
iments on the dynamics of inhomogeneous turbulence in 
4 He (vibrating grid) [3] and 3 He-B (propagating turbu- 
lent front) [15| hinted at new dissipation mechanisms in 
T = limit too. 

The need for experiments with homogeneous turbu- 
lence in superfluid 4 He generated by large-scale flow in 
T = limit has led us to an experiment in which energy- 
containing eddies comparable in size to that of the con- 
tainer are generated by an impulsive spin-down to rest of 
a rotating cell of a square cross-section. It is known that 
with classical fluids spin-down to rest is always unstable 
resulting in 3-dimensional turbulence after a few radians 
of the initial rotation The values of L in our ex- 

periments, monitored by firing a beam of charged probes 
of known trapping diameter through the tangle, were in 
the range 10 1 — 10 5 cm" 2 , corresponding to characteristic 
inter-vortex distances £ of 3 - 0.03 mm. 
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FIG. 1: (color online) Cartoon of the vortex configurations 
in the experimental cell (side view) at different stages, (a) 
Regular array of vortex lines during rotation at constant Q 
before deceleration, (b) Immediately after stopping rotation 
(0 < t < 10/Q,), turbulence appears at the outer edges but 
not on the axis of rotation, (c) After about t ~ 30/il, 3D 
homogeneous turbulence is everywhere, (d) At t ~ 10 3 /fi 
most of the 3D turbulence has decayed. Shaded areas indicate 
the paths of probe ions when sampling the vortex density in 
the transverse (L t , c, red) and axial (L a , d, green) directions. 



The inner experimental volume of the container had 
the shape of a cube of 4.5 cm side; it was filled with pure 
4 He at a pressure of 0.1 bar, and mounted on a rotating 
cryostat (Fig. 1; for details see [13] )• Negative ions could 
be injected by either of the two tungsten field-emission 
tips [18| through gridded holes in the centers of the bot- 
tom and of one of the side plates. Then the ions could 
be pulled through the volume by an applied electric field 
E (converging to counteract the space-charge repulsion) 
and collected by the electrodes in the center of the op- 
posite plates protected by Frisch grids. To make sure 
the free decay of the tangle and measurements of L are 
not affected by the space charge, only one short pulse 
(duration < 0.5 s) of injected ions was fired, at time t af- 



ter stopping rotation. Thus each data point in Figs. 2-4 
represents a different realization of the turbulence. Be- 
fore decelerating to rest, the cryostat was kept in steady 
rotation at the required Vt for at least 300 s. The decel- 
eration was linear in time, taking 2.5 s for Q = 1.5 rad/s 
and 0.1 s for f2 = 0.05 rad/s. The origin t = was chosen 
at the start of deceleration. 
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FIG. 2: (color online) L t (t) at T = 0.15 K for four values of 
fi. Average electric fields used for Q = 1.5 rad/s: 5 V/cm 
(o), 10 V/cm (A), 20 V/cm (o), 25 V/cm (v). The dashed 
line shows the dependence t~ 3 ^ 2 . Horizontal bars indicate the 
equilibrium values of L at = 1.5, 0.5, 0.15. 0.05 rad/s (from 
top to bottom). 



The pulse of collector current arrives after a well- 
defined time of flight across the cell, and the rela- 
tive reduction in its amplitude, 7(t)/J(oo), is converted 
into the average vortex density through L(t)/Lo = 
(crcQ -1 ln(J(oo)/J(i)). There are two types of charge car- 
rier |l|. Free ions, dominating at T > 0.8 K, are electron- 
containing bubbles of radius « 19 A; at T — 1.6 K their 
trapping diameter a oc E^ 1 was 0.1 /im at E = 20 V/cm, 
and the time of flight was 1 s. The others, dominating at 
T < 0.7 K, are small quantized vortex rings of diameter 
D ~ 1 /im with one electron bubble trapped in its core; at 
T = 0.15 K their trapping diameter er ~ D was between 
0.4 and 1.7 /xm for E between 5 and 25 V/cm. In the nar- 
row range of temperatures , 0.7-0. 8 K , both types can ex- 
ist, distinguishable at sufficiently high fields because the 
free ions move much faster than the rings; here measure- 
ments of L with both types of carriers agreed well. Their 
trapping diameters a and times of flight were calibrated 
in situ, for all temperatures and field strengths used, on 
arrays of parallel rectilinear vortex lines with a range of 
equilibrium densities L = 2Q/k produced by steady rota- 
tion at different an gula r velocities f2. They agreed with 
the published data [19|, [201 • We could not measure val- 
ues of L that are too high (no detectable current reaches 
the collector) or too low (no detectable change in cur- 
rent). For charged vortex rings and free ions, the accessi- 
ble ranges of L were about 10 1 -10 4 and 5 • 10 2 -10 6 cm" 2 , 
respectively. At low temperatures, with Kelvin waves ex- 



cited across a range of wave numbers, the vortex length 
is a fractal property, and hence can be underestimated 
by the use of probe objects of size a ~ 1 fim. However, 
as the values of L at T — 0.15 K, measured for a range 
of a = 0.4-1.7 /im, were in good agreement, any unde- 
tected contribution to L must be small. In Fig. 2, the 
measured densities of vortex lines along the horizontal 
axis (transverse, L t ) are shown for four different initial 
angular velocities O. During the transient, which lasts 
~ 10/ O, L t (t) goes through a maximum, after which it 
decays, eventially reaching the universal late-time form 
of L oc t- 3 / 2 . 
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FIG. 3: (color online) fT 3/2 L t (t) (filled symbols) and 
S7 -3 / 2 L a (£) (open symbols) vs. D.t for four values of Q at 
T = 0.15 K. The straight line oc t~ 3/2 guides the eye. Hori- 
zontal bars indicate the equilibrium values of L at = 0.05, 
0.15, 0.5. 1.5 rad/s (from top to bottom). 



In Fig. 3 the measured densities of vortex lines along 
the horizontal, L t , and vertical (axial, L a ) axes are 
shown, by solid and open symbols respectively. To stress 
the scaling of the characteristic times with the initial 
turn-over time O -1 and the universal late-time decay 
oc t~ 3 / 2 , the data for different O are rescaled accord- 
ingly. At all O the transients are universal. Immediately 
after deceleration, L a is stable near the equilibrium den- 
sity at rotation, 20/ K, while L t jumps above this value, 
indicating the appearance of turbulence at the perimeter. 
Only at t rs 3/0 does the former start to grow, signalling 
the destruction of the rotating core consisting of vertical 
rectilinear vortices. After passing through a maximum 
at t = 8/0 and t — 15/0 respectively, L t and L a merge 
at t ~ 30/O and then become indistinguishible. This im- 
plies that from then on the turbulence in the cell becomes 
isotropic and probably homogeneous. The scaling of the 
transient times with the turn-over time O -1 tells us that 
transient flows are similar at different initial velocities O, 
as expected for flow instabilities in classical liquids. 

At 0.08 < T < 0.5 K the observed L(t) and hence in- 
ferred v'{T) were independent of temperature. In Fig. 4, 
we compare the transients L t (t) at low (0.15 K) and high 



1<T 



6 10 J 



a 

^ ION 



10' 



T- 1 {.V^** 

1 = 1 .6 K % 



0.15 K 



* 0.05 rad/s 
0.15 rad/s 
0.5 rad/s 

° 1.5 rad/s 




10" 



10" 



10' 



io z 



10' 



10* 



at 



FIG. 4: (color online) Q,- 3/2 L t {t) vs. Qt for T = 0.15 K 
(filled symbols) and T — 1.6 K (open symbols). Dashed and 
solid lines oc t -3 / 2 guide the eye at T = 1.6 K and 0.15 K, 
respectively. 



(1.6 K) temperatures. Their durations (positions of the 
maximum in L{t)) are virtually identical. On the other 
hand, at T — 0.15 K, the late-time decay L oc i~ 3 / 2 
takes longer to develop (time ~ 100/O) and then the 
prefactor is about 10 times larger than for T = 1.6 K. 
This implies that at low temperatures the steady-state 
inertial cascade with a constant energy-containing eddy 
size and constant energy flux down the range of length 
scales requires a much greater total vortex line density 
and perhaps some extra time to build up this density. 
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FIG. 5: (color online) The effective kinematic viscosity v' 
after a spin-down from Q = 1.5 rad/s measured in the trans- 
verse (A) and axial (v) directions. Closed (open) triangles 
correspond to measurements with free ions (charged vortex 
rings). Error bars specify the uncertainty of fitting. Squares 
and diamonds: second sound measurements of grid turbulence 
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Our values of v\ along with those obtained with towed 
grids of two different designs in a channel of 1 x 1 cm 2 
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cross-section at high temperatures (corrected by the fac- 
tor of 0.72 as suggested in [22j]), are shown in Fig. 5. In 
discussing the behaviour of v'(T), we note first that in 
the low temperature limit a quasi-classical inertial-range 
energy spectrum cannot persist to wave numbers greater 
in order of magnitude than n/£. At larger wave numbers 
quantum effects must dominate. Therefore the maxi- 
mum vorticity associated with the quasi-classical spec- 
trum cannot be greater than that given by 

<W 2 >^x S) « / k 2 E k dk = ^1^/3^/3^2/3 (4) 

Jo 4 
The energy flux e is given by Eq. JI]), therefore 

Now k 2 L 2 can be identified as an effective total mean 
square vorticity in the system Q. Thus, if u'/k < 0.09 
the quasi-classical spectrum cannot account for all the 
vorticity in the system. In the low-temperature limit 
some, but not, according to our estimates, all, of the ex- 
tra vorticity could be associated with the Kelvin-wave 
cascade; some must exist either at the high-wavenumber 
end of the quasi-classical cascade, associated with a bot- 
tleneck, or in an intermediate regime of the type de- 
scribed by Kozik and Svistunov 0]. According to the 
detailed calculations of L'vov et al. Q the ratio v'/n is 
only about (ln(£/ao))~ 5 = 10 -6 , and therefore our mea- 
surements suggest that any bottleneck is less severe than 



these calculations predict. We remark that, although 
there is vorticity in excess of that given by the Kol- 
mogorov cascade, the associated energy must be small 
enough not to invalidate our assumption that the largest 
eddies account for most of the energy; otherwise the i~ 3 / 2 
decay would not be observed. 

The striking temperature dependence oiv'/n shown in 
Fig. 5 therefore suggests to us the following scenario. As 
the temperature falls from 1.3 K to 0.8 K the Kolmogorov 
cascade expands to reach wavenumbers of order tt/£. Be- 
low 0.8 K (when the estimated mutual friction parameter 
a ~ p n /p becomes a < 10 -3 ) the flow of energy for a 
given line density drops sharply as vorticity accumulates 
at wavenumbers around tt/£ and Kelvin waves begin to 
be excited. As the temperature falls further the damping 
of Kelvin waves by mutual friction decreases until below 
0.5 K the Kelvin-wave cascade extends to wavenumbers 
far beyond ir/£. The low value of v' in the zero tem- 
perature limit is a measure of the difficulty in transfer- 
ring energy through wavenumbers around ir/£ from the 
3-dimensional quasi-classical Kolmogorov cascade to the 
1-dimensional Kelvin-wave cascade. Better understand- 
ing of this process remains a challenge for theorists. 
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